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Abstract 

Optimization procedures allow one to design a 
spur gear reduction for maximum life and other end 
use criteria. A modified feasible directions 
search algorithm permits a wide variety of inequal- 
ity constraints and exact design requirements to be 
met with low sensitivity to initial guess values. 
The optimization algorithm is described and the 
models for gear life and performance are presented. 
The algorithm is compact and has been programmed 
for execution on a desk top computer. Two examples 
are presented to illustrate the method and its 
application. 

Nomenclature 

A velocity factor 

B material surface constant, psi 

b design parameter scaling constant 

C design problem constant 

C dynamic capacity, lb 

d design parameter scaling coefficient 

E elastic modulus, psi 

F force, lb 

/ face width, in. 

Vf feasible direction gradient vector 

Vh total violated constraint gradient vector 

J AGMA bending stress geometry factor 

6 

f life (10 cycles or hours) 

M merit function 

Vm merit function gradient vector 

N number of teeth 

R gear radius, in. 

R reliability 

RMS surface finish 

S ac allowable surface strength, psi 

As optimizing step size 

T temperature, 0 F 

V constraint 

V velocity, in. /sec 

V stress volume, in. 3 

Vv constraint gradient vector 

W weight, lb 

X unsealed design vector 


X AGMA factor 

Y scaled design vector 

z depth to critical stress, in. 

a angle of approach, deg 

6 characteristic life ( 10 6 cycles or hours) 

p friction coefficient 

V Poisson's ratio 

p surface radius of curvature, in. 

g contact stress, psi 

T q shear stress, psi 

<p pressure angle, deg 

W angular velocity, rad/sec 

Exponents 

b Weibull slope 

c proportionality factor 

h proportionality factor 

p load-life exponent 

Subscripts 

10 90 percent probability of survival 

( reliability) 

ag gear addendum 

B base 

b bending 

c center distance 

d total dynamic 

f flash 

g gear 

H Hertzian 

i independent parameter index 

j search step index 

k constraint index 

Li parameter lower estimate 

pi pinion radius to gear addendum on line of 
action 

s sliding 

t tooth 

Ui parameter upper estimate 

w weight 



Introduction 


The optimal design of a spur gear mesh is a 
problem of considerable interest in mechanical 
design. 1 ’ 4 When fracture of the gear teeth due to 
bending is the primary mode of failure, the minimum 
number of teeth which avoids interference offers 
the strongest gear set for a given size. However, 
as speeds increase, so do the prospects for pitting 
and scoring modes of failure. Pitting at and below 
the pitch point on the pinion tooth is a failure 
mode which limits the life of the gear teeth. 
Procedures have been presented to design a gear set 
for minimum size and minimum weight, considering ^ 
pitting fatigue as well as bending and scoring. 
These procedures determine an optimal gear set from 
a weight or size standpoint alone. 

One promise of optimization is adaptability to 
the application at hand. For example, the shaft 
center distance may be fixed by other considera- 
tions in the design of a machine. Given this situ- 
ation, one would like to apply optimization theory 
to determine a design which maximizes the gear set 
life for a given center distance. 


as specified by upper and lower bounds on the inde 
pendent design parameters, X, such that: 


- 1.0 < Y A < 1.0 


( 1 ) 


as 

x li < X i < X ui 

Thus : 

■ d ! X l * b l 

where 
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and 


(5) 


The computer is useful here, and the accessi- 
bility of the personal computer makes this use 
attractive. A modified feasible directions search 
algorithm is applied in a continuous design space 
in a form which is memory and time efficient. 

Using a fixed search step, which is halved only 
when the search passes a solution, simplifies the 
computations. The method is an improvement of the 
basic gradient search method. 7,8 It checks the 
gradients in the inequality constraints as well as 
the merit function to calculate a feasible search 
direction which improves the merit function while 
staying within the acceptable design region. After 
finding the continuous optimum design, the program 
allows the designer to check alternate designs 
which may be more practical. These designs can 
have parameter values which obey the discrete 
parameter limitations and are close to the contin- 
uous optimal design's parameter values. 

This method is applied to the problem of 
designing for maximum Life at a given center dis- 
tance as an example. Also, it is applied to a sim- 
ilar problem of minimum center distance with a 
given life requirement for comparison purposes. 


Optimization Method 

As with most optimization techniques, the pro- 
cedure begins with several vectors. These vectors 
are the independent design variables, X? the 
inequality constraints, V; the parameters of the 
merit function, P; and the constants which define 
the specific problem, C. An optimization solution 
is the design variable values, X, which minimize or 
maximize the merit function value while maintaining 
all constraint values, V, inside their specified 
limits. A procedure starts with a guess for the 
design variable, X, and iterates with some logical 
procedure to find the optimal design variable. 

Scaling 

To maintain balance among the independent 
design parameters, the design space is scaled into 
a continuous, dimensionless design space. The 
scaled design parameters, Y, vary from -1.0 to 1.0 


The actual design variable, X^ can always be 
retrieved from the scaled variable, Y^, by: 



Gradient 

Central to the method is the gradient calcula- 
tion. This is performed with small perturbations 
in the design variables from the nominal position. 
The gradient in the merit function, Vm, is calcu- 
lated as: 

0M 

0M 

3 ^ 


0y 3 


where, 

0M 

M<Y J ,...,Y 1 .-AY,...,*,,) - M<Y 1( ... ,Yj, Y„) 

S* (8) 

In Eq. (8), M is the merit function, Y l is a 
scaled design variable, and A Y is the small 
change which is made in each Y . In the program, 
Ay is set at 0.001 which is 0.05 percent of the 
full range of a scaled design parameter. 
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The magnitude of the gradient vector is given 

by i 


I Vm | 


E 

i-i 


5m 


(9) 


For minimization, the direction of change in Y 
which reduces the merit function, M, at the great- 
est rate is determined by the unit vector, Vm: 


VM 



( 10 ) 


For maximization, the sign in Eq. (10) reverses. 

In the simple gradient method, Eq. (10) 
defines the direction for the step change in the 
scaled design vector. 

Y j + 1 - Yj + As Vm (11) 


in the constraint value, V^, for motion in the 
gradient direction. The vector sum of the gradi- 
ents in the violated constraints, Vh, is the second 
gradient of the algorithm: 


EK 

Vh - _!! 

£K 


(14) 


The gradient in the violated constraints, Vh, 
point 3 towards the acceptable design space from the 
unacceptable design space. By itself, it enables 
the algorithm to turn an unacceptable initial guess 
into an acceptable initial guess by a succession of 
steps : 

Y j+1 - Y^ + As Vh (15) 


where As is the scalar magnitude of the step. If 
no constraints are violated, this will be the next 
value for Y in the search. 

Step Size 

Step size, As, is a significant element of any 
optimization procedure. 8 For stability and direct- 
ness, a fixed step size is used. Initially, the 
step size is 0.1, which is 5 percent of the range 
of a single design parameter. But the procedure 
halves the step whenever a local minimum is reached 
or the search is trapped in a constraint corner. 

To complete the simple procedure, the search 
declares a solution when the percent change in the 
merit function, M, is less than a pre-set limit of 
0 , 0001 . 


|M J*i - M ;i 


< 0.0001 


( 12 ) 


Feasible Direction 


Once inside the acceptable design region, the 
algorithm proceeds along the steepest descent 
direction until the calculated step places the next 
trial outside the acceptable design space. To 
avoid this condition, the algorithm selects a fea- 
sible direction for the next step. Figure 1 shows 
a sloped constraint intersecting contour lines of 
improving merit function values. The figure shows 
gradients in the merit function, Vm, and the 
impending constraint, Vh. The two gradient vectors 
are placed at the last viable design step, although 
the constraint gradient is calculated at the trial 
location - one step in the merit gradient direction 
from the indicated point. The feasible direction 
selected, Vf, is the unit vector sum of these two 
gradients in the merit function, Vm, and the vio- 
lated constraints, Vh: 


Vm + Vh 
f?m + Vh| 


(16) 


Initial Value 

The optimization procedure described above is 
scaled, fixed step, and steepest decent. When the 
initial guess is in the acceptable design space, 
and the optimum is a relative minimum, this method 
works quite well. However, placing the initial 
guess in the acceptable design region is often a 
problem, and in many cases, the best design is 
determined by a "trade-off" among conflicting 
design constraints at the edge of the design space. 

These problems are addressed with a second 
gradient. Just as one can calculate the gradient 
in the merit function, one can calculate the gradi- 
ent in a constraint variable: 



where Vv is a unit vector in the direction of 
decreasing value in the constraint, V^. For upper 
bound constraints, moving through the design space 
in the direction of Vv^ will reduce the con- 
straint value V . For lower bound constraints, a 
sign reversal in Eq. (13) will produce an increase 


And the next step becomes: 

Yj 4l - Y 3 + As Vf ( 17 ) 


Algorithm Use 

By using subroutines to calculate the merit 
function and constraint values for each design 
trial, the procedure separates the logic of the 
algorithm from the analysis necessary to define the 
problem. This allows the design problem to be 
changed easily without concern for the optimization 
procedure. The directness of the procedure adds 
additional steps, but enables the program to run on 
a personal computer. 

An additional benefit of separating the anal- 
ysis routines from the optimization logic is the 
ability to modify the design at execution and ver- 
ify the characteristics of similar, more practical 
designs with the same program. The optimization 
procedure works with a continuous design space, 
which includes gears with fractions of teeth and 
nonstandard sizes. By allowing the user to see the 
ideal continuous variable solution and to modify 
this to designs with whole numbers of teeth and 
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standard sizes, the procedure enables a designer to 
determine a practical optimum design easily. 

Input Data 

The algorithm requires the user to provide an 
ASCII data file with four groups of data: the con- 

stants, design variables, constraints, and merit 
function weighting coefficients. All constants, 
variables, limits and coefficients have names and 
units in the data frle to make the program output 
more understandable. In addition, the data file 
includes a title for the optimization task. 

The constants define the specific problem out- 
side of the analysis subroutines. With their 
labeling names and units, these constants define 
the specific problem being optimized in the output 
data file. The constants also enable similar 
problems to be solved by merely changing the con- 
stant values without recompiling and linking the 
analysis subroutines to the optimizing routines. 

The design variables are the independent 
parameters for which the algorithm seeks values. 
Three values are included for each variable: a 

lower estimate, an upper estimate, and an initial 
guess. The lower and upper estimates determine the 
scaling range as the span between them is set to 
two dimensionless units by the optimizing proce- 
dure. However, they do not set hard limits for the 
search. If hard limits are needed to avoid 
indeterminate calculations, these may be set in the 
constraint limits. 


file at each step to document the search path and 
assure the user that it is running properly. 

Once an optimum is found, the program writes 
the design variable values that produced the opti- 
mum, the merit function value and its component 
parameter array, and the full constraint variable 
array to the screen and the output file. It pauses 
for user inspection, and on resumption gives the 
user a chance to enter a different set of design 
variable values. This allows the user to include 
practical, near optimal design trials in the com- 
puter analysis record. Since the best design 
should be of this type, this provides full documen- 
tation to the additional trials. When the user 
declines to enter more trials, the program closes 
the output file and ends execution. 


Analysis Subroutines 

Two analysis subroutines apply the optimi- 
zation procedure to the 9pur gear design problem. 
The routines are subroutine BOUNDS and subroutine 
VALUES. Subroutine BOUNDS analyzes a gear design 
for its constraint variable values. Subroutine 
VALUES analyzes a gear design for its merit func- 
tion parameter values. These subroutines and an 
input data file treat the two problems of designing 
for maximum life at a given center distance and of 
designing for minimum center distance at a given 
life. 

Bounds 


The constraint limits include a bound value 
and a direction: lower or upper. Both a lower and 

an upper bound may be set on the same property by 
including two separate limits in the constraint 
list. Since the program reports the constraint 
values for the optimal and check designs, one can 
add inactive constraints which are always satisfid 
to obtain a printout of additional properties in 
the constraint list. This also gives one the 
opportunity to switch the controlling constraints 
to obtain different designs for different physical 
conditions without changing the program. 

The last group of data are the parameters 
which may be included in the merit function and 
their weighting coefficients. These weighting 
coefficients may convert different properties to a 
common measure such as cost. They may be order of 
magnitude corrections to keep one parameter from 
dominating the design. Or they may all be zero 
except one, to select that parameter as the 
quantity to be optimized, so the same program can 
be used to obtain different designs for different 
requirements by switching the merit function 
parameter with the non-zero coefficient. This 
switch changes the objective of optimization, while 
the switch of active constraints changes the 
environment of optimization. 

Runtime Options 

In running the program, the user is requested 
to enter the data file name prefix. The program 
will look for this name on an input data file with 
the suffix ".IN" and will write the output to a 
file with this prefix and the suffix ".OUT" as well 
as to the screen as it runs. Since the program 
calls the analysis subroutines many times, it 
writes information to the screen and the output 


Included in the 16 constraint variables evalu- 
ated by subroutine BOUNDS are: some distances, the 

pinion weight and loading, and some stress, life 
and scoring limits. The involute interference var- 
iable is the distance along the line of action from 
the base circle of the pinion to the addendum cir- 
cle of the gear. This distance must be positive 
for the gear tooth tip to contact the pinion tooth 
on its involute surface and avoid interference. 


The pinion weight is the product of the 
density of steel times the volume of a disk with a 
diameter equal to the pinion pitch diameter and a 
thickness equal to the gear and pinion face width. 
The transmitted load is the pinion torque divided 
by the base radius of the pinion, which is the 
nominal force acting between the gears along the 
line of action. The pitch line velocity is the 
rotational speed of the pinion times its pitch 
radius . 


In this analysis, the dynamic load estimate is 
the AG HA velocity factor model. In terms of a gear 
quality number, Q , the AGMA estimate of the sum of 
the transmitted load and the dynamic load is: 


F 


d 


A + /v 
\ A , 


(18) 


where 


A * 50 + 56 1 


(12 - Q v ] 


(19) 
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In Eq. (18), F is the total dynamic load, F t is 
the nominal transmitted load and V is the pitch 
line velocity of the gears. In Eq. (19), the gear 
quality number, Q , may have a value between 6 and 
11 with 11 corresponding to the higher quality 
gear. All gear stresses and lives are calculated 
using this total dynamic load, with a quality 
number, Q * 9. 

V 


the load sharing factor, X M is a thermal - elastic 
factor and X Q is a geometry factor. 

The pinion life and mesh life calculations are 
described in the gear life model section which 
follows . 

Values 


Gear tooth bending fatigue and gear tip scor- 
ing are modes of failure considered in the program. 
The bending fatigue model uses the AGMA J factor to 
estimate the bending stress with the load at the 
highest point of single tooth loading on the pin- 
ion. The formula for the bending stress is: 

F • P 

0 , If (20) 

b / ‘ J 

where J is the AGMA J factor 14 and / is the 
effective gear tooth face width. 

The maximum contact stress and gear tip 
Hertzian pressure are calculated as: 



l/ f? 

♦ x/ h | 

7T f COS ^ 

2 

1 - *V 




+ : 


E 

E„ 


L p 

9 )\ 


The. second analysis subroutine, VALUES, calcu- 
lates the three components of the merit function: 
the mesh life in thousand hours, the center dis- 
tance in inches, and the pinion weight in pounds. 
The total merit function is: 

M - c,l + C c (R p + R g ) + C w w (24) 

where 9 is the mesh life, the radii sum is the 
center distance and W is the pinion weight. The 
three weighting function coefficients: and 

C are constants defined in the input data file. 

By assigning different values to these coeffi- 
cients, the user can change the problem being 
solved. With not equal to zero and the other 

two coefficients equal to zero, the optimization 
will seek out a solution which maximizes the gear 
mesh life. With and C w equal to zero and 

C negative, the optimization will seek out a 
solution which minimizes the center distance. 


Gear Life Model 


where the p * s are the radii of curvature of the 
two gear tooth surfaces at the point of contact, 
the V s are the Poisson ratios for the two gear's 
materials and the E's are the elastic moduli for 
the two gear's materials. The maximum contact 
stress occurs at the lowest point of full load con- 
tact on the pinion tooth. The gear tip Hertzian 
pressure uses one half of the total dynamic load 
since the load is shared between two tooth pairs at 
this point. 

The gear tip scoring model includes the pres- 
sure times velocity factor and the critical oil 
scoring temperature model from lubrication theory. 
The normal pressure times sliding velocity is pro- 
portional to the frictional power loss of the gear 
set. This factor is the highest for contact at the 
gear tip, so the normal pressure is the gear tip 
Hertzian pressure. the sliding velocity at the 
gear tip is given by: 

V . - W g R .g 3in It + a .9> " W pV 3in + <V (22) 


The gear life model, based on surface pitting, 
comes from rolling element bearings. 6 Surface pit- 
ting of gear teeth follow a similar pattern to 
bearing race pitting, with the possible difference 
of surface initiation. Lundberg and Palmgren pro- 
posed the model in the late 1930 's. 10 They assumed 
that the log of the reciprocal of the reliability, 
R, of a bearing is proportional to its life, 9 , and 
some stress parameters. These parameters are: the 

stress level, T ; the depth to the maximum shear 
stress, z ; and the stress volume, V. The rela- 
tionship is: 



T z 


-h 

O 


V| 


(25) 


In relation (25) b is the Weibull slope and 
c and h are exponents of proportionality to be 
found experimentally. This is the equation for the 
two-parameter Weibull distribution with the addi- 
tion of stress and size factors. In general terms, 
the two-parameter Weibull distribution is: 


where the W's are the angular velocities of the 
gears, the R's are the radii to the contact point 
on the two gears, $ is the nominal pressure angle 
of the gear set and the a's are the angles of 
approach on the two gears. 

The lubricating oil flash temperature is 
another factor used to monitor gear tooth scoring. 
One estimate of this temperature 14 is given by: 


LN 


(26) 


where b is the Weibull slope or shape factor and 
9 is the characteristic life of the distribution. 
To replace the characteristic life with a 90 per- 
cent probability of survival life, f iQ , solve 
Eq. (26) for the characteristic life, 6. 
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0.45 
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( 27 ) 


where T 0 is the base temperature of the oil, jl a 
is the average surface friction coefficient, X^ is 
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At R » 0.9, the life is l iQ . Substituting 
this into Eq. (26) gives the two parameter Weibull 

life as: 


(28) 


in terms of 

the 
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/ \ 
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The life to reliability relationship of 
Eq. (28) is for a specific load which determines 
the 1 life. This load, F, is related to the 
component dynamic capacity, C, ass 


* io 



(29) 


The dynamic capacity of the whole gear is 
lower than that of a single tooth. In a single 
pass, fixed axis gear set, each rotation of the 
gear subjects every tooth on the gear to a single 
load cycle. The gear fails when any single tooth 
on the gear fails, and the fatigue damage in each 
tooth accumulates independently of the damage in 
the other teeth. In successive coin tosses, the 
probability of a specific combined event is the 
product of the probabilities of each coin toss. So 
too, the reliability of the gear, R , is the 
product of the reliabilities of each tooth in the 
gear . 


R 


g 


(33) 


Here, the dynamic capacity of the component, C, is 
the load which has a 90 percent reliability life of 
one million cycles, the load on the component is 
F, and the power, p, is the load-life exponent. 
Since the life at the dynamic capacity is one 
million load cycles or unity, it does not appear as 
a variable in the equation. 

Because gear tooth life behaves similar to 
bearing life, engineers at the NASA Lewis Research 
Center formulated a model for gear tooth life simi- 
lar to the bearing life model. 6 ' 11 Starting with 
Eq. (25) Coy, Townsend and Zaretsky developed a 
model for the reliability and life of a spur gear. 
The model uses both the two-parameter Weibull dis- 
tribution of Eq. (28) and the Palmgren load-life 
relation of Eq. (29), With statistically repli- 
cated data, they showed that these models predict 
gear tooth pitting. 

From Eq. (25), they determined a relationship 
for the dynamic capacity, C^, of a spur gear tooth. 
Rounding the exponents to whole numbers gives the 
dynamic capacity as a function of Buckingham's 

. . , . _ 12,13 

load-stress factor, B. 


c t 


\u/p) 


(30) 


In Eq. (30), B is a material strength which has 
the dimensional units of stress, the effective face 
width of the tooth is /, and the curvature sum at 
the failure point for the contacting teeth is 
El//?. The curvature sum is: 

ll/p - 1 + -L (31) 

p 9 p p 


Here, p is the radius of curvature of the 
gear tooth surface at the failure point, and p 
is the radius of curvature of the pinion tooth 
surface at the failure point. With the dynamic 
capacity expressed in this form, the material 
strength factor serves the role of the surface 
fatigue strength, S , of the AGMA design code. 

A relation for the material strength factor in 
terms of the surface fatigue strength is: 



In Eq. (33), N is the number of teeth on the 
gear, and R 9 is the reliability of a single tooth 
on the gear. The reliability of any tooth on a 
gear is equal to the reliability of any other tooth 
on the gear. To transform Eq. (33) into a life 
relationship, substitute Eq. (26) for the two reli- 
abilities into the log of Eq. (33). 


I 


10, g 



10, t 


(34) 


The gear life, i , has units of million 
io , g 

gear rotations. 


The gear dynamic capacity, C , is found by 
substituting Eq. (34) into Eq. (2§) for the tooth 
to obtain the analog of Eq. (29) for the gear. 
This produces: 


C 


9 


N 


l/(b*p) 

9 


(35) 


Gear Design 

Consider the design of a gear set to transmit 
80 hp from a shaft turning at 5000 rpm to an output 
shaft turning at 2500 rpm. The life of the gears 
as a pair in hours should be maximized at a center 
distance of 5.0 in- 

Table 1 is the program output listing of the 
problem defining constants and the design variable 
range and initial value. The gears are to have a 
20® nominal pressure angle and be made of high 
strength, heat treated steel with a tooth surface 
finish of 32 rms. The material surface constant of 
9800 psi corresponds to a surface compression 
endurance strength of 200 000 psi at 10 7 fatigue 
cycles and a reliability of 90 percent. The load- 
life factor of 8.93 is from the ANSI/AGMA 2001 B88 
Standard, 14 and the Weibull slope of 2.5 is from 
the NASA Lewis gear test data. 6/11 The three 
design variables to be found are the number of 
teeth on the pinion, the diametral pitch and the 
face width of the gears. For the gear, the number 
of teeth will be twice that of the pinion. 

The gears should have a face width which gives 
the pinion a length to pitch diameter ratio between 
0,2 and 0.5. Additional strength limits placed on 
the design are: a tooth bending stress limit of 
40 000 psi, a surface contact stress limit of 
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150 000 psi, a pressure times velocity factor of 
100 million psi-ft/min, and an oil flash temper- 
ature limit of 275 °F. Involute interference 
should be avoided as well. These limits are shown 
in Table 2 which lists the design constraints and 
merit function weighting factors in the output data 
file. 

The controlling constraints are the request 
for a center distance no larger than 5.0 in. and 
the maximum pinion length to diameter ratio. 
Included in the list of constraints are properties 
which have limits that any design will satisfy, 
such as: pinion torque, transmitted load, dynamic 

load, pinion life, and the oil flash temperature. 
These constraints enable the program to analyze 
each design for their values and display them as 
part of the design summary for each optimum and 
check design. A velocity factor equation calcu- 
lates the dynamic load as prescribed by ANSl/AGMA 
2001 B88 with a quality number of 9. All calcu- 
lated stresses and contact pressures use the 
dynamic load. 

Table 3 summarizes the program's optimal solu- 
tion for the maximum life design including the 
design variable values, the merit function param- 
eter values and all constraint values and limits. 
For thi9 ideal design, the mesh life is 28 000 hr. 
Table 4 gives a similar summary for a selected 
maximum life design. Figure 2 shows this design 
which has 46 and 92 teeth on the two gears, a 
diametral pitch of 14 and a face width of 1.625 in. 
It has a life of 18 000 hr with all other con- 
straints satisfied. The design has a pinion weight 
of 3.9 lb, a transmitted load of 650 lb and a total 
dynamic load of 1840 lb. 

To illustrate the flexibility of the design 
program, the design is changed to request a minimum 
life of 2000 hr to see how compact the design can 
become. 

Six changes are made to the input file. The 
center distance limit changes from an upper bound 
of 5.0 in. to a lower bound of 0.0 in. The mesh 
life limit changes from a lower bound of 0.0 hr to 
a lower bound of 2000 hr. The merit function 
direction changes from MAX to MIN. And, the non- 
zero weighting factor switches from mesh life to 
center distance. Labeling changes for the file 
itself and the problem title make up the last two 
changes which enable the program to preserve and 
distinguish the input and output records for the 
two design requests. 

Table 5 summarizes the ideal optimum design 
for the minimum size design. In this design, the 
gear center distance is reduced to below 5.5 in. by 
the request to reduce the mesh life by an order of 
magnitude. Table 6 summarizes the nearest practi- 
cal design, which is shown in Fig. 3. In this 
design, the pinion and the gear have 36 and 72 
teeth, and both have a diametral pitch of 12 and a 
face width of 1.5 in. The center distance is 
4.5 in. and the mesh life is 2200 hr. 

Although the bending and contact stresses have 
increased over those for the maximum life design, 
they are still within the acceptable criteria. The 
scoring limit values also have increased to accept- 
ably larger values. Had any of these constraints 
been reached, the bounding constraints would have 
influenced the ideal and final optimal designs. 


It should be noted that all the designs have 
the maximum possible face width. Designs based on 
pinion weight rather than center distance also 
required the maximum possible face width. Minimiz- 
ing pinion weight for a life of 2000 hr produced 
the same optimal design as that found by minimizing 
the gear center distance. 

Conclusions 

By combining the power of optimization with 
the access of the desk top computer, a practical 
spur gear life design program has been written. 

The program uses a fixed step, modified feasible 
directions algorithm to search for the optimum 
design. Basic relations for the optimization 
algorithm are presented. 

Extensive labels and keyboard interactions 
give the designer a record of the ideal optimal 
design and any user specified designs entered after 
the optimum has been found and reported. Small 
changes in the input data file enable one to redi- 
rect the design objective from maximum life to 
minimum size and to change the controlling design 
constraints . 

Two gear analysis subroutines and the compat- 
ible input data file apply the optimization proce- 
dure to the gear life design problem. The first 
subroutine analyzes a gear design for the con- 
straint variable values. And the second subroutine 
analyzes a gear design for the merit function 
parameter values. The input data file defines 
these values, sets the limits for the design con- 
straints, and selects the active merit function 
parameters . 

Tooth surface pitting fatigue life produces 
the finite life of the gear set. A fatigue life 
model for this mode of failure is presented. Sta- 
tistical variations in gear life as predicted by 
the two parameter Weibull distribution are given 
also. 

Two designs illustrate the design calculations 
of the program. One gear set is designed for maxi- 
mum life with a required center distance, and a 
second design is obtained for a minimum center 
distance at a minimum acceptable life. The design 
procedure considers other constraints as well to 
obtain practical gear designs. Minimum weight and 
minimum center distance designs were found to be 
identical . 
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TABLE 1. - GEAR DESIGN CONSTANTS AND VARIABLES 
FOR SPUR GEAR LIFE MAXIMIZATION 

[Design with modified gradient optimization using a maximum step 
limit and scaled variables.] 

(a) Fixed design requirements 


Poisson's ratio 0.25 

Elastic modulus, psi 3x10 

Pressure angle, deg 20 

Gear ratio ...... 2 

Transmitted power, hp ....... 80 

Pinion speed, rpm 5000 

Material weight density, lb/in. 3 0.283 

Material surface constant, psi , 9800 

Weibull slope 2.5 

Load-life factor 8.93 

Reliability ........... ... 0.9 

Base temperature, °F ..... ...... 120 

Tooth surface fiish, rms 32 


(b) Estimated values of the three independent 
design variables 



Low 

High 

Initial 

Pinion teeth 

10 

100 

40 

Diametral pitch, in. -1 

4 

28 

14 

Face width, in. 

0.5 

5 

2.5 
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TABLE 2. 


16 MAXIMUM LIFE CONSTRAINTS AND MERIT FUNCTION 


(a) The 16 constraint limits 


Constraint 

Limit 

Type of 
bound 

Involute interference, in. 

0.001 

Lower 

Lower face width-to-diameter ratio 

0.2 

Lower 

Upper face width-to-diameter ratio 

0.5 

Upper 

Pinion weight, lb 

0 

Lower 

Center distance, in. 

o 

in 

Upper 

Pinion torque, lb-in. 

0 

Lower 

Transmitted load, lb 

0 

Lower 

Total dynamic load, lb 

0 

Lower 

AGMA bending stress, psi 

0.4xl0 5 

Upper 

Full load contact stress, psi 

1.5xl0 5 

Upper 

Gear tip hertz pressure, psi 

1.5xl0 5 

Upper 

Pinion life, cycles 

0 

Lower 

Mesh life, hr 

0 

Lower 

Pitch line velocity, ft/min 

0 

Lower 

PV factor, psi-ft/min 

100 

Upper 

Flash temperature, °F 

275 

Upper 


(b) Maximize the objective function (OBJ), which 


has three terms. OBJ * the linear sum of 


Component 

Unit 

Multiplied 

by 

Mesh life 

in thousands of hours 

1 

Center distance 

in inches 

0 

Pinion weight 

in pounds 

0 


TABLE 3. - MAXIMUM LIFE OPTIMAL DESIGN 
(Optimization successful in 50 steps. ] 


(a) The final design vectors 


Design parameter 

Scaled 

vector, 

Actual 

vector, 

Pinion teeth 
Diametral pitch, in." 
Face width, in. 

-0.24857 

-.23797 

-.48148 

43.81443 

13.14433 

1.66667 


(b) Components of the maximum objective 
function, 28.1361 


Component 

Value 

Multiplied 

by 

Mesh life, hrxlO 3 

28.136 

1 

Center distance, in. 

5.0000 

0 

Pinion weight, lb 

4.1161 

0 


(c) The 16 constraint values 


Constant 

Value 

Limit 

Type of 
bound 

Involute interference, in. 

0.23224 

lxlO’ 3 

Lower 

Lower face width-to-diameter ratio 

0.5 

0.2 

Lower 

Upper face width-to-diameter ratio 

0.5 

0.5 

Upper 

Pinion weight, lb 

4.1161 

0 

Lower 

Center distance, in. 

5.0 

o 

in 

Upper 

Pinion torque, lb-in. 

1008.4 

0 

Lower 

Transmitted load, lb 

643.87 

0 

Lower 

Total dynamic load, lb 

1825.9 

0 

Lower 

AGMA bending stress, psi 

0.34422x10 s 

0.4xl0 5 

Upper 

Full load contact stress, psi 

1.2290xl0 S 

1.5x10 s 

Upper 

Gear tip hertz pressure, psi 

0.80590x10 s 

1.5x10 s 

Upper 

Pinion life, cycles 

9. 5274xl0 9 

0 

Lower 

Mesh life, hr 

28 . 136xl0 3 

0 

Lower 

Pitch line velocity, ft/min 

4100.2 

0 

Lower 

PV factor, psi-ft/min 

67.268xl0 6 

100x10 s 

Upper 

Flash temperature, °F 

203.96 

275 

Upper 
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TABLE 4. - MAXIMUM LIFE DESIGN CHECK 


(a) Design check 


Design parameter 

Scaled vector, 

x i 

Pinion teeth 

46.000 

Diametral pitch, in. 

14.000 

Face width, in. 

1.625 

(c) The 


(b) The components of the maximum objective 
function, 18.1249 


Component 

Value 

Multiplied 

by 

Mesh life, hrxlO 3 

18.125 

1 

Center distance, in. 

4.9286 

0 

Pinion weight, lb 

3.8993 

0 


values 


Constraint 

Value 

Limit 

Type of 
bound 

Involute interference, in. 

0.23810 

lxlO" 3 

Lower 

Lower face width-to-diamater ratio 

0.49457 

0.2 

Lower 

Upper face width-to-diameter ratio 

0.49457 

0.5 

Upper 

Pinion weight, lb 

3.8993 

0 

Lower 

Center distance, in. 

4.9286 

5.0 

Upper 

Pinion torque, ib-in. 

1008.4 

0 

Lower 

Transmitted load, lb 

653.20 

0 

Lower 

Total dynamic load, lb 

1843.0 

0 

Lower 

AGMA bending stress, psi 

0.37648xl0 5 

0 . 4xl0 5 

Upper 

Full load contact stress, psi 

1.2582x10 s 

1,5x10 s 

Upper 

Gear tip hertz pressure, psi 

0.81064x10 s 

1.5xl0 5 

Upper 

Pinion life, cycles 

6 . 1374xl0 9 

0 

Lower 

Mesh life, hr 

18 . 125xl0 3 

0 

Lower 

Pitch line velocity, ft/min 

4041.6 

0 

Lower 

PV t factor, psi-ft/min 

64. 271xl0 6 

100x10® 

Upper 

Flash temperature, °F 

202.37 

275 

Upper 


TABLE 5. - MINIMUM SIZE OPTIMAL DESIGN 
[Optimization successful in 27 steps. ] 


(a) The final design vectors 


Design parameter 

Scaled vector. 

Actual vector, 
x < 

Pinion teeth 
Diametral pitch, in." 1 
Face width, in. 

-0.45045 

-.36590 

-.55778 

34.72992 

11.60917 

1.49501 


(b) The components of the minimum objective 
function, 4.48739 


Component 

Value 

Multiplied 

by 

Mesh life, hrxlO 3 

2.0236 

0 

Center distance, in. 

4.4874 

1 

Pinion weight, lb 

2.9739 

0 


(c) The 16 constraint values 


Constraint 

Value 

Limit 

Type of 
bound 

Involute interference, in. 

0.20055 

lxlO" 3 

Lower 

Lower face width-to-diameter ratio 

0.49919 

0.2 

Lower 

Upper face width-to-diameter ratio 

0. 49919 

0.5 

Upper 

Pinion weight, lb 

2.9699 

0 

Lower 

Center distance, in. 

4.4870 

0 

Lower 

Pinion torque, lb-in. 

1008.4 

0 

Lower 

Transmitted load, lb 

717.48 

0 

Lower 

Total dynamic load, lb 

1965.3 

0 

Lower 

AGMA bending stress, psi 

0.38256xl0 S 

0 . 4xl0 5 

Upper 

Full load contact stress, psi 

1.4325x10 s 

1.5x10 s 

Upper 

Gear tip hertz pressure, psi 

0.98539x10 s 

1.5x10 s 

Upper 

Pinion life, cycles 

0 . 67728xl0 9 

0 

Lower 

Mesh life, hr 

2 . OOOlxlO 3 

2 . 0xl0 3 

Lower 

Pitch line velocity, ft/min 

3679.5 

0 

Lower 

PV factor, psi-ft/min 

92 . 242xl0 6 

lOOxlO 6 

Upper 

Flash temperature, °F 

234.28 

275 

Upper 
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TABLE 6. - MINIMUM SIZE DESIGN CHECK 


(a) The final design vectors 


Design parameter 

Actual, vector, 

X . 

Pinion teeth 

36.000 

Diametral pitch, in.' 

12.000 

Face width, in. 

1.500 


function, 4.500 


Component 

Value 

Multiplied 

by 

Mesh life, hrxlO 3 

2.1964 

0 

Center distance, in. 

4.5000 

1 

Pinion weight, lb 

3.0006 

0 


(c) The 16 constraint values 


Constraint 

Value 

Limit 

Type of 
bound 

Involute interference, in. 

0.20588 

lxlO" 3 

Lower 

Lower face width-to-diameter ratio 

0.5 

0.2 

Lower 

Upper face width-to-diameter ratio 

0.5 

0.5 

Upper 

Pinion weight, lb 

3.0006 

0 

Lower 

Center distance, in. 

4.5 

0 

Lower 

Pinion torque, lb-in. 

1008.4 

0 

Lower 

Transmitted load, lb 

715.41 

0 

Lower 

Total dynamic load, lb 

1961.4 

0 

Lower 

AT3MA bending stress, psi 

0 . 38980xl0 5 

0.4x10 

Upper 

Full load contact stress, psi 

1 . 4239xl0 5 

1 . 5x10* 

Upper 

Gear tip hertz pressure, psi 

0.97093xl0 5 

1 . 5xl0 5 

Upper 

Pinion life, cycles 

0.74373xl0 9 

0 

Lower 

Mesh life, hr 

2 . 1964xl0 3 

2 . 0xl0 3 

Lower 

Pitch line velocity, ft/min 

3690.2 

0 

Lower 

PV factor, psi-ft/min 

88.071x10® 

100x10 s 

Upper 

Flash temperature, °F 

230.41 

275 

Upper 



Figure 1 . — Gradient sum to find feasible search direction. 
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